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Abstract 

We investigate an undamped random phase-space dynamics in deterministic external force 
fields (conservative and magnetic ones). By employing the hydrodynamical formalism for those 
stochastic processes we analyze microscopic kinetic- type "collision invariants" and their relation- 
ship to local conservation laws (moment equations) in the fully nonequlibrium context. We address 
an issue of the continual heat absorption (particles "energization") in the course of the process 
and its possible physical implementations. 

PACS number(s): 02.50.-r, 05.40.-a, 03.65.-w 

1 Introduction 

The major topic discussed in the present paper will be the frictionless random motion which is in- 
troduced in close affinity with second-order processes driven by white noise, 012- I n the course of 
motion an unrestricted "energization" of particles is possible, a phenomenon which has not received 
much attention in the literature, although appears to be set on convincing phenomenological grounds 
in specific (nontypical) physical surroundings, [310]. 

A classic problem of an irreversible behaviour of a particle embedded in the large encompassing 
system, jSUHJEl pertains to the random evolution of tracer particles in a fluid (like e.g. those performing 
Brownian motion). In that case, a particle - bath coupling is expected to imply quick relaxation 
towards a stationary probability distribution (equilibrium Maxwellian for a given a priori temperature 
T of the surrounding fluid/bath) and thus making the system amenable to standard fluctuation- 
dissipation theorems. The dissipation of randomly acquired energy back to the bath is here enforced by 
dynamical friction mechanisms, 151151 111)1 171. although no balance between stochastic forcing (random 
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acceleration effects) and dissipative losses of energy is in fact in existence in the course of the process, 
unless asymptotically, [ffllllj. 

The traditional picture of a tracer particle in a fluid/gas, even if modeled in terms of stochastic 
diffusion-type processes, is usually considered in direct reference to the Boltzmann kinetic theory. That 
justifies an exploitation of properly defined microscopic collision rules to yield dynamical friction (c.f. 
0^]) and then to allow for a "derivation" (which breaks down the microscopic conservation laws 
on the way, of the Brownian dynamics in suitable scaling limits. The crucial assumption about 
a rapid decay of correlation functions amounts to strong friction and short relaxation time, 0E1IE1- 

Let us mention that, as reveals the hydrodynamical formalism of the Brownian motion, |14| . 
Brownian particles "while being thermalized" , show up a continual net absorption of heat/energy 
from the reservoir (that point is normally disregarded in view of the very short relaxation times): so- 
called kinetic temperature grows up to the actual temperature of the bath, in parallel with the mean 
kinetic energy of the Brownian particle, ^2 ED E3- c.f. also [121 E|. Indeed, J7| chap. 2, for t > 
we encounter on the average a continual " heating" phenomenon in the course of which an asymptotic 
equilibrium value kT/m is achieved: i?|l?| 2 = %[1 — exp(— 2/3t)] — ► ^ where q — f3kT/m, while it 
stands for velocity random variable. In this case the bath is thought of as insensitive to this process 
and is supposed to remain in the perpetual state of "statistical rest". 

(Notice that such situation is surely in conflict with dynamical approaches to the thermalization 
of the Brownian particle, ^0E|j which take into account that the particle-thermostat coupling acts 
both ways (consider reaction forces, third Newton law in the mean etc.). In consequence, it is both 
the Brownian particle and the dissipating agent (bath, thermostat) which need to thermalize together 
i.e. simultaneously evolve towards a common state of statistical equilibrium, cf. 

In below we shall consider purely classical prototype stochastic models which depart from most 
traditional Brownian motion in its Langevin (white noise), Kramers and Smoluchowski versions. How- 
ever our focus will be on an extremally nonequilibrium situation when no thermalisation is possible. 
Then, stationary (equilibrium) state is definitely out of the reach, no fluctuation-dissipation relation- 
ships are in existence and the untamed heating phenomenon is the most conspicuous characteristic of 
the randomly forced dynamics. 

1.1 Phenomenological motivations 

In contrast to the well established frictional " Brownian motion standards" , we shall explore the 
opposite regime of dissipation - free theory. The departure from the standard theory appears to be 
quite radical, since it refers to a different physics, where dissipative time scales (if any) are much 
longer than the time duration of processes of interest, including the particle life-time. 

We are particularly motivated by two papers E] devoted to the thermal energization of particles 
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which is due to impulsive stochastic forcing. In fact, frictionless stochastic processes were invoked to 
analyze situations present in magnetospheric environments. Specifically, one deals there with charged 
particles in a locally uniform magnetic field which experience inetrmittent stochastic electrical forcing. 

However, then a serious conceptual obstacle is present, if the theory is considered without suitable 
reservations: one needs to be able to control and eventually tame the kinetic energy growth, since in 
the absence of friction we deal with an efficient particle "energization" mechanism. That involves as 
well an important issue (not addressed in the present paper) of suitable thermalization/thermostatting 
mechanisms that would possibly tame the otherwise unlimited "energization" of particles or make it 
effectively irrelevant on the time-scales of interest. 

Small random impulses occur within ambient fields in such way that the smooth deterministic 
dynamics of a particle is irregularly disrupted by sudden and very short bursts of low-intensity noise. 
Without friction, we deal with an undamped random walk and no fluctuation-dissipation relationship 
can be established. 

Phenomena of that kind are seldom mentioned in the traditional statistical physics research, since 
they may be realized in rather unusual, from the point of view of the common practice, environments. 

Indeed, |2()j the physical nature of environments in which the random propagation is investigated, 
appears to be crucial. In the original papers E] , the environment of interest was an astrophysical or 
space plasma, where densities and therefore dissipation are often reduced by many orders of magnitude. 
That amounts to dissipation time scales which are much longer than those on which any random forcing 
may take place. Moreover, physical processes associated with the thcrmalization were eliminated in 
view of the fact that particles entered geometrical regions of space where the plasma instabilities 
(producing stochastic electric fields fluctuations) and thus the source of randomness had disappeared. 
In contrast to typical laboratory environments (e.g. plasmas), the high spatial heterogeneity present 
in space plasmas makes random forcing to occur only within very limited regions of the non-stationary 
environment and not to occur everywhere and for all times. 

That means that it is the very nature of the environment which causes the source of energization 
to disappear and reappear spontaneously, while typical random bath models heavily rely on the noise 
spatial homogeneity and its "eternal" existence, [Sj. 

1.2 Hydrodynamical formalism for random dynamics 

To elucidate the formal roots of our strategy and reasons for studying the frictionless random dynamics 
by means of the hy drodynamical formalism, let us call our attention back to the standard Ornstein- 
Uhlenbeck (dissipative) framework. The Langevin equation for mass m particle in an external field of 
force (we consider a conservative case) F = F (~x) = — W reads: 
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q. = -^+£+2(0 (2 ) 

Random acceleration A (t) obeys the white noise statistics: (Ai (s)) = and (Ai (s) Aj (s 1 )) = 
2qS (s — s') Sij, where i = 1,2, 3. 

Since things are specialized to the standard Brownian motion, we know a priori that noise intensity 
is determined by a parameter q — D0 2 where D — , while the friction parameter is given by the 
Stokes formula m[3 = 6nr]a (or its analogue in case of the Lorentz gas, [xj). Consequently, the effect 
of the surrounding medium on the motion of the particle is described by two parameters: friction 
constant (3 and bath temperature T. Assumptions about the asymptotic (equilibrium) Maxwell- 
Boltzmann distribution and the fluid reaction upon the moving particle are here implicit, 

The resulting (Markov) phase - space diffusion process is completely determined by the transition 
probability density P ( ~x, it, t\ ~xo, ~Tio, to), which is typically expected to be a fundamental solution 
of the Kramers equation: 



9P n 
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V^P - q^P (3) 



The associated spatial Smoluchowski diffusion process with a forward drift b (~x) = is analyzed 
in terms of increments of the normalized Wiener process W(t). The infinitesimal increment of the 
configuration (position) random variable X(t) reads: 

dX (t) = ^—dt + V2DdW (t) (4) 

mp 

The related Fokker-Planck equation for the spatial probability density p(~x, t) reads dtp = DAp — 
V • (pb) and explicitly employs the large friction regime, |10llUlH3*]. In fact, we take for granted that 
both time and space scales of interest (i.e. those upon which the accumulation of relevant random 
events prove to be significant) largely exceed the relaxation time interval /3 _1 and that dominant 
contributions "of interest" come from velocities |lx | < (kT/m) = {q/P) 1 ^ 2 and that the corresponding 
variation of "r* is sufficiently small (actually it is of the order I'ul/P = (<?//3 3 ) 1/ ' 2 ), Pl- 
under those assumptions the Fokker -Planck equation for the spatial Markov process arises as 
the scaling (/? > 1) limit of the 0-th order moment equation associated with the original Kramers 
law of random phase-space dynamics. In fact, by following the traditional pattern of hydrodynamical 
formalism, (141 116| . we infer the closed system of two (which is special to Markovian diffusions !) local 
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conservation laws for the Smoluchowski process, [111118) : 



d tP +V-(l?p) = (5) 
(d t +-u-W)~v = V(fi-Q). (6) 

Here (we use a short-hand notation ~Tf(~x,t) = ~u) 

^0?,i) = ^--D^ (7) 
mp p 

defines so-called current velocity of Brownian particles and, when inserted to Eq. (5), transforms the 
continuity equation into the Fokker- Planck equation, |10) . 

Eq. (6) stands for the scaling limit of the first order moment equation derivable from the kinetic 
equation (3) and directly corresponds to the familiar Euler equation, characterizing the momentum 
conservation law in the lowest order approximation of kinetic theory based on the Boltzmann equation, 

ng. 

However, the large friction regime enforces here a marked difference in the local momentum con- 
servation law, in comparison with the standard Euler equation for a nonviscous fluid or gas. Namely, 
instead of the kinetic theory motivated expression for e.g. rarified gas: 

. -> ~F VP 

(d t + l?-V)l?= (8) 

m p 

where P(l?) stands for the pressure function (to be fixed by a suitable equation of state) and F is the 
very same (conservative - VF) force acting upon particles as that appearing in the Kramers equation 
(3), the Smoluchowski regime (6) employs the emergent volume force (notice the positive sign) +VS1 
instead of — V^: 

2\mPJ ymfi J w 

and the pressure-type contribution — VQ where, j^I] (see also [211 \12\ ) 

Q = 2D 2 ^JT- (10) 

and A = V 2 is the Laplace operator, does not leave any room for additional constraints upon the 
system (like e.g. the familiar equation of state). 

To have a glimpse of a dramatic difference between physical messages conveyed respectively by 
equations (6) and (8), it is enough to insert in (8) the standard equation of state P(x) — ap 13 with 
a,f3 > and choose F = —oj 2 ~x to represent the harmonic attraction in Eqs. (2) - (9), see also [T3] , 
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Markovian diffusion processes with the inverted sign of V (CI — Q) in the local momentum conser- 
vation law (6) i. e. respecting 



(9 t +u*-V)l? = V(Q-fi) (ff) 

instead of Eq. (6), were considered in Ref. as implementations of the "third Newton law in 
the mean". Nonetheless, also under those premises, the volume force term — \7Q in Eq. (11) does 
not in general coincide with the externally acting conservative force contribution (e.g. acceleration) 
±F = -J-W akin to Eq. (8). 

Accordingly, the effects of external force fields acting upon particles are significantly distorted 
while passing to the local conservation laws in the large friction (Smoluchowski) regime. 

That becomes even more conspicuous in case of the Brownian motion of a charged particle in the 
constant magnetic field, ^HHIEH]- In the Smoluchowski (large friction) regime, friction completely 
smoothes out any rotational (due to the Lorentz force) features of the process. In the corresponding 
local momentum conservation law there is no volume force contribution at all and merely the " pressure- 
type" potential Q appears in a rescaled form, [T5] : 

a 2 Afl 1 / 2 
Q = WT^f 2D ^ (12) 

where (3 is the (large) friction parameter and ui c = is the cyclotron frequency of the charge q e 
particle in a constant homogeneous magnetic field B = (0,0, B). Clearly for moderate frequency 
values uj c (hence the magnetic field intensity) and sufficiently large (3 even this minor scaling remnant 
of the original Lorentz force would effectively disappear, yielding Eq. (10). 

This observation is to be compared with results of Refs. 0E] where, in the absence of friction, the 
rotational Lorentz force input clearly survives when passing to the local conservation laws, in plain 
contrast with the Smoluchowski regime. By disregarding friction (or considering it to be irrelevant 
on suitable time scales) it is possible to reproduce exactly the conservative external force acting upon 
particles in the local conservation laws, |24) . 

Thus the low noise regime, if combined with the sufficiently short duration time of processes of 
interest (much lower than the dissipation time scales) sets the range of validity of local conservation 
laws in the classic Eulcr form (8), when demanded to comply with the stochastic diffusion-type 
microscopic dynamics assumption. 



2 Random motion of a free particle 
2.1 Resume 



Let us consider the Langevin-type equation for a particle which is not suffering any friction while 
being subject to random accelerations 



dx 
dt 

du ... 



(13) 



We fix initial conditions: x(t ) = Xq, u(to) — Uq, t a = 0. The fluctuating term A (t) is as- 
sumed to display a standard white noise statistic: (i) A (t) is independent of u, (ii) (A (s)) = and 
(A (s) A (s 1 )) = 2qS (s — s'), where q is here left as an unspecified (albeit a priori physical) parameter. 

This random dynamics induces a transition probability density P (x,u, t\ xq, uq, to) which uniquely 
defines the corresponding phase-space Markovian diffusion process executed by (x,u). Here, the 
function P (x, u, t\ xq, uq, to) is the fundamental solution of the Fokker-Planck equation associated 
with the above Langevin equation: 



' u — + 1l^2- ( 14 ) 



dP _ dP d 2 P 
dt dx ^ du 2 

A fundamental solution of Eq. (14) was first given by Kolmogorov and is a literature classic, 
For the uniformity of further argumentation, we shall here follow a direct solution method 
[261 IT7| which employs properties of the involved Wiener process Wt ■ 

Since \/2qWt — f A (s) ds, the distribution of the random variable S = u — uq is gaussian with 
mean zero and variance 2qt. We define the displacement of the particle as x — xo — f u (rj) drj. Hence, 
we have R — x — xo — uot = \[2q f Q W s ds. 

Because S and R are gaussian random variables, their joint probability distribution W (R, S) is also 



gaussian and is given in terms of mean values and covariances of S and R. The mean values are equal 

to zero, the covariance matrix is given that C = ( C ° V ft % C ° V ft g )=2 9 (K 

\ cov(b,R) cov(b,b) J \ ±t 

cov{W Sl W t ) = s At. 

Thus,M/(i?,5) = 5 L( a ^)^exp r 1 ' " ' ' ' 1 U 



S 



\{R S )C- 1 

Since the problem of finding the transition probability distribution P (x,u, t\ xq, uq, to) is equivalent 
to finding the joint probability distribution W (R, S) of random variables S and R, we have thus arrived 
at the following expression for (time homogeneous) P (x,u, t\ Xo, uo, to)' 



P(x,u,t\ x ,uo,t = 0) = —— 2 exp 



1 VT2 
2^2qT 2 



(u-uo) 2 3(x-x ~^t) 21 



(15) 



Aqt qt 3 

We encounter a situation when standard fluctuation-dissipation relationships are manifestly invalid 
Nonetheless, it is useful to mention a useful computational connection with the familiar damped case 
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Namely, let us consider the Langevin equation for the free (would be Brownian) particle with the 
friction term 

I- < 16 » 

^ = -(3u + A(t) (17) 

but no a priori presumed (fluctuation-dissipation) relationship between (3 and q. We again fix the 
initial conditions: x (to) — xq, u (to) — uq, to = 0. 

Except for considering (3 and q as independent parameters, we can directly exploit the classic result 
[2] for the corresponding transition probability density P (x,u,t\xo,uo,to = 0) solving the Kramers 
equation (compare e.g. Eq. (3)) 

dP dP n d , . d 2 P 

-m +u dx- = ^ Pu) + q l^ (18) 

Its solution is known to arise in the form of a joint probability density W (i?, S) of random variables 
S = u - uoe' 131 and R = x — x — /3 _1 (l — e _/3t ) u Q 



gR 2 - 2hRS + fS 2 



2(fg-h 2 ) 

where the auxiliary coefficients are given that f = jjs (2/3t — e -2 "* + 4e - "* — 3j, g = j| (l — e _2/3t ) 
h=j r (l-e-»)>. 

Since we assume that the intensity of noise q is independent of the friction parameter (3, we can 
safely take the limit of (3 — > in the joint probability density W (i?, S 1 ). Then, in view of the obvious 
limiting values for S = v — vo, R — x — xo — Vot, f — f^t 3 , <? = 2gt, /i = qt 2 we readily arrive at the 
previous frictionlcss-case expression. 

Let us stress that it is the assumption about the independence of parameters q and f3 which takes 
us away from the standard theory of Brownian motion. We cannot any longer expect that thermal 
equilibrium conditions can be approached by standard dissipation mechanisms. We also recall that 
in case of the "normal" dissipative free Brownian theory its physical formulation is based on the 
requirement that the distribution of velocities must be given by the stationary Maxwell-Boltzmann 
probability density with variance as t — > oo. Moreover, the traditional Einstein-Smoluchowski 
theory relates the spatial diffusion coefficient D to (3 by the fluctuation - dissipation relation D = 
That appears to be consistent with the Kramers picture of the phase-space random dynamics only if 
we a priori assume that q=^(3 1 [HI El- Then q = D(3 2 and large (3 definitely implies large q. 

2.2 Kramers equation and local conservation laws 

The stationary (time homogeneous) transition probability density for the diffusion process without 
friction in phase-space has the form (15) and is the fundamental solution to Kramers (Fokker- Planck 



(19) 
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type) equation (14). We are interested in passing to a hydrodynamical picture, following the traditional 
recipes ^] E] ■ To that end we need to propagate certain initial probability density and investigate 
effects of the random dynamics. Let us choose most obvious example of: 

nM = 4 « P (-<£^>!) (J^)% IP (_<^) . (20) 

so that at time t we have p (x, u,t) — J P(x,u, t\ xq,uq, to = 0) po (xq, uq) cLxqCLuq. 

Since P (x, u,t\ xq, uq, to) is the fundamental solution of the Kramers equation, the joint density 
p(x, u,t) is also the solution and can be written in the familiar, j^j, form of Eq. (19) for p(x, u, t) = 
W(R,S). However, in the present case functional entries are adopted to the frictionless motion and 
read as follows: S = u — Uj n j, R = x — Xi n i — Ui n it with 

/ = a 2 + b 2 t 2 + ^qt 3 

g = b 2 + 2qt (21) 
h = b 2 t + qt 2 . 

The marginals p(x,t) — J p (x, u, t) du and p (u, t) = J p (x, u, t) dx are 



and 

^ V > \2ng) V \ 2gJ \2k (b 2 + 2qt) ) l \ 2 (b 2 + 2qt) ) V ' 

Let us introduce an auxiliary (reduced) distribution W {S\R) = j w(SR)ds wnere m the denomi- 
nator we recognize the marginal spatial distribution J W (S, R)dS = w of W(S, R). 

Following the standard hydrodynamical picture method [141 lltjl 1 1 1 j we define local (configuration 
space conditioned) moments: (u) x = J uWdu and (u 2 \ x — J u 2 Wdu. First, we compute 

h b 2 t + qt 2 

j a z + bH A + =:qt 6 



h 2 \ _ qt 3 (2b 2 + qt) + 3 a 2 (b 2 + 2qt) 



3a 2 + < 2 (3b 2 + 2qt) (25) 
Now, the first two moment equations for the Kramers equation are easily derivable. Namely, the 
continuity (0-th moment) and the momentum conservation (first moment) equations come out in the 
form 
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dw d ,. . . 

^«u>,t») + ^«u a >.»«;) = 
These equations yield the local momentum conservation law in the familiar form 
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d 



1 ftP*. 



u> dx 

where we encounter the standard |16| textbook notion of the pressure function 

Pkin{x,t)= (u 2 ) x ~(u)l w(x,t) . 



The marginal density w obeys = — 2/V 



Aw 1 



1 dPkin 

w dx 



= 2 (fg - h 2 ) V 



and that in turn implies 
Aw 1 / 2 ' 



.1/2 



As a consequence, the local conservation law Ij28(l takes the form 



^ + <u> x V)(uL = 



= +2(/ 5 -^ 2 )V 



Aw 1 / 2 



,,1/2 



+VQ 



where (we point out the plus sign in the above, see e.g. Eq. (11)) 



(26) 
(27) 

(28) 

(29) 



(30) 



(31) 



fg-h 2 = a 2 b 2 + 2a 2 qt + h 2 qt 3 + ^q 2 t* = D 2 (t) 



(32) 



and by adopting the notation D 2 (t) = fg — h 2 we get — ^ = +VQ with the functional form of 
Q given by Eq. 10. Here, instead of a diffusion constant D we insert the (positive) time -dependent 
function D(t). 

With those notational adjustments, we recognize in Eq. (31) a consistent Euler form of the local 
momentum conservation law, in case of vanishing volume forces (c.f. Eq. (8) and compare that with 
Eq. (11) again). 

It is useful to mention that by employing Eq.(23) we can readily evaluate the mean kinetic energy 
associated with the frictionless free dynamics. Namely, there holds 



< u(t) 2 > - < u(t) > 2 = b 2 + 2qt 



(33) 



where < u(t) >—< u(0) >= u ini . That is to be compared with the standard Brownian motion 
expression, ^7j: < u(t) 2 >= %[1 — exp(— 2/3t)]. In the present case we definitely encounter an 
untamed, continual heating (energization of Refs. 0E1) mechanism. 

Let us observe that limits a — > 0, b — > are here under control and allow to reproduce a consistent 
hydrodynamical formalism for the transition probability density of the frictionless process. Indeed, in 
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that case, we readily reproduce observations of Refs. that < u 2 (t) >= 2qt, while for the spatial 

mean displacement we would get < x 2 (t) >— §<Z^ 3 - The spatial process shows features typical of an 
anomalous transport and obeys the Fokker-Planck equation dtp = (qt 2 )Ap. 

3 Charged particle in a constant magnetic field 

Our analysis of a free particle subject to random accelerations and no friction can be extended to the 
case when external deterministic force fields are in action. As a specific example we shall discuss the 
dynamics of a charged particle in a constant magnetic field, see e.g. also |15l 1231 |2~§1 I3(J| . 

Magnetic field effect upon the charged particle in practically frictionless random environments 
was the main objective of Refs. 0] in the formulation of a kinematic mechanism for particle 
energization. Those authors were primarily interested in the application of their model in concrete 
physical situations, hence they did not concentrate on finding the general solution to the involved 
probabilistic problem. They have evaluated the mean square velocity and gave computer assisted 
visualization of the dynamics. However no general descriptions of the related velocity space nor 
phase-space stochastic process were given. Thus no information was available about spatial features 
of the dynamics, nor about the associated macroscopic (hydrodynamical formalism) balance equations. 

3.1 Velocity-space process 

For a charged particle in a constant magnetic field and fluctuating electric environment, we consider 
the equation of motion in the form 

^ = il-tf X B + A(t) (34) 
at mc 

where it denotes the velocity of the particle of an electric charge q e and mass m. A fluctuating 
part A (t) represents random forces of electric origin which are to obey the standard white noise 
statistics assumptions (c.f. Section 2.1). 

Let us assume for simplicity that magnetic field B is directed along the z-axis of a Cartesian 
reference frame so that B = (0, 0, B) and B = const. In this case, Eq. (|34(1 takes the form 

iif. = -Ku + A it) (35) 
at 

where 

/ -uj c \ 
A= I lu c 
\ / 

and ui c = denotes the cyclotron frequency. 
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The corresponding transition probability density P ( u , t\ u o, to = Q) of the time-homogeneous 
Markov process with the initial condition: P (~u, t\ u o) — * S 3 (u — ~uq) as t — > has the form: 

P(^^o)=(^j ex P (-J (36) 

where U(t) = exp(— At), see e.g. Ref. [T5] . 

P (~u , i|l?0i ^0 = 0) is the solution of the Fokker-Planck equation: 

dP 

— = -u e [V^ x PH]^ + qV^P (37) 

where [V^» x P"u] i=3 = {Pu2) — gf^ (Pui). The expression (|36|l can be readily extended to the 
case t = s ^ 0, t > s, 

P ^^ S ^(^(^)j eXP (~ 4 g (U) J" (38) 
Let us consider in some detail the motion of an individual particle which has an initial velocity ~uq 
at tg = 0, and the initial probability density g ( v , s = 0) = go (1?) = 5 (v — ~??o)- The propagation 
in time is controlled by the transition density distribution P (it , s) according to 

Q(Tt,t) = j ' P(~u,t\~v,s = ti)Q {~v)dv (39) 

8 &>*> = U*J eXP ( y - J i^J' (40) 

By evaluating mean values of the resultant density for i = 1, 2 velocity components, we obtain: 



(tti(t)) = / iiig(u ,t)du (41) 

(l?(t)) = (K(t)),(u 2 (t))) = C/(t)l?o (42) 

/oo 

u 2 iQ {!t,t)dlt (43) 
-OC 

: <^ 2 W)-( Ml (i)) 2 = 2< Z t (44) 



<7, 



we realize that the kinetic energy of the particle increases linearly with time: (u 2 {t)) — (u 2 (0)) = Aqt. 
This result obviously coincides with less explicit calculations of 0. 

If we define the displacement of the particle as 7* - ~ro — L ~u (rj) drj, we can extend the 
description of the process to the phase space. An attempt to specify the transition probability density 
P (~r^, ~u , ilV'o, ~uo, to — 0) is equivalent to finding the joint distribution W (^S, R^j of random vectors 

— > — > 

S and R , which are defined by the equations 

~S = ~u-U(t)lt (45) 
~R = -"r -A _1 [l-C/(t)j u (46) 
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where A 



° ^ I and C/(t) 





The joint probability density of S and R has the form 



cos oj c t sin w c i 
- sin uj c t cos w c t 



IT 



1 



A-K 2 {fg-h 2 -k 2 ) 
with the auxiliary coefficients: 



exp 







2 






/ 


~s 


+ g 


R 


2 - 2hS -~R +2k (l) x ~R^j 









2(fg-h 2 - k 2 ) 



g = 2qt 



(47) 



t sin u! c t 

(Jr. 



2q , 

h = — ^ (1 — cosuj c t) 
(Ji 



(48) 



2q 



1 . 

t sin ui c t 

U>r 



Eq. il4"?|) is the solution to the Fokker-Planck equation in the phase space: 



dW 



(49) 



where [V^ x WTt]^ = ^ (Wu 2 ) - £ (Wm). 
3.2 Kramers equation and local conservation laws 

The formula i|47|) is an explicit solution to the Kramers Fokker-Planck equation and coincides with 
the phase-space transition probability density P ( r , ~u , t\ ~r^o, uq, to = 0). 
Let us consider the following initial phase-space probability density 



1 



2ttw 



• exp 



2a 2 



1 



2vr& 2 



exp 



(1?0 - ~UiniY 

2b 2 



(50) 



which is propagated in the course of the stochastic process so that at time t we have p(r , u , i) = 
/ P ( "r , "u , i| 7*0, ^o, k = 0) Po C^0) ~Uo) dV (Tu 

Since P ( r , u ,t\ r q, Uo,to = 0) is the fundamental solution of the Kramers equation the joint 
density p ( r , u , t) is also the solution. This joint density can be put into the form 



p(?,lt,t) = w(~S,l?) 



(51) 



4tt 2 (fg -h 2 - k 2 ) 



exp 



S 



R 



2hS ■ R+2k( S x R 



i=3 



2{fg-h 2 - k 2 ) 
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where 



S = ~U —U(t) U ini 

R = ~r - ~r lnl - A" 1 [1-U (*)] ~u im 

g = b 2 + 2qt 

2b 2 u c (1 — cos u c t) + a 2 uj^ + 4g (u> c t — smiv c t) 

h 

k 



2q ( 1 — cos uj c t ) + b 2 uj c sin w c t 



2q (ui c t — sinuj c t) + b 2 uj c (1 — cos oj c t) 



The marginals p (V, t) = J p ( r , u , t) d~u and p (~u , t) = J p ( r , w , t) rfV read 

2^ 



p(r ' <) = 2^/ exp 



1? 



2/ 



(52) 
(53) 
(54) 
(55) 

(56) 

(57) 



(58) 



2n-^j (2b 2 uj c (1 — coso; c t) + a 2 u>^ + Aq {u) c t — sinw c t)) 



x exp 



- r im - a- 1 [i-u (t)\ 



and 



2-K; (2b 2 u> c (1 — cos Lo c t) + a 2 uj3 -f 4q (uj c t — sinu> c t)) 
I 



p(^,t) = o — exp 
2irg 



s 



\ 



2.9 



2ir (b 2 + 2qt) 



exp 



I u -U (t) u ini [ 
2 (b 2 + 2qt) 



(59) 



Notice that in view of < u 2 (t) > — < Ui(t) > 2 = b 2 + 2qt, where < Ui{t) >= [U(t)~Ui n i\i with 
i = 1,2, and < ~u(t) > 2 = \U{t)~Ui n i\ 2 — \~Ui n i\ 2 —< ~u(0) >, we have the same as before (see e.g. 
Section 2) indicative of the linear in time energy growth: 



< ~u 2 (t) > - < ~u(0) > 2 = 2b 2 + 4qt . 



(60) 



We shall perform major steps of the hydrodynamical analysis following the pattern of Sec. 2. To 

w( !>,~r) 



this end we introduce the auxiliary distribution W 



/_» _>/ __> which has an explicit form 

/ wl S ,R\dS 



W 



P^) = 271 



l(fg-h 2 -k 2 ) 



exp 



S — m 



V 



2±{fg-h 2 -k 2 ) 



(61) 



where m = j (hRi — kR 2 , hR 2 + fci?i), R and S are given by equations (52) and (53) respectively. 

We define local moments as (ui)-^ ~ J UiWcTu and (uf)-^ = J ujWdu where i = 1,2. It is 
apparent from the Eq. (UH that (~u)-fi = (( u i)-fi > ("2)^) = U (t) ~Ui n i + m and (uf)-^ = (ui) 2 ^ + 

j(fg-h 2 -k 2 ). 



14 



The first two moment equations for the Kramers equation (|49(l are easily derivable: 



d t w + V [(u)^ w] = 

d d 

dt[(ui)^w]+ — [(u 1 )^w]+ — [(u 1 )^(u 2 )^w] = u c (u 2 )^w 



Q 

dt [{u 2 )^ w] + ^ [(u 2 )-g (u 2 ) 3 H 

Those equations imply the local momentum conservation law 



u)-a = -A (It) V-*P 



— ft?)-* x B--V ■ p" 



(62) 
(63) 

(64) 



(65) 



— > ■* — * . . . < dP < — * 

where V ■ P is a vector which i-th component is equal to 2_^. = 1>2 and P denotes the 

tensor consisting of components Py = ((ui ~ ( u i)^J ( u j ~ they are given that Pn = 
((u?>3 - (ux)|) u>, P 22 = ((n|>-g - (u 2 )|) w and P 12 = P 21 = 0. 



Let us point out that 



(66) 



displays a linear growth in time ~ qt for large i. 
Owing to relation 

ldPu idP 22 



w dx ' w dy 
the local momentum conservation law l|65|l takes the form 



2 (/ 5 -h 2 - k 2 ) V 



Aw 1 / 2 



d t + (It)^ V 



0^ x B +2(fg~h 2 -k 2 ) V 



Aw 1 / 2 



where the crucial time-dependent coefficient reads 
fg — h — k = a b \- + f - + -2-=- - 



2a qiH j- cos(tuj c ) — sm(tui c 



(67) 



(68) 



(69) 



In contrast to the free dynamics coefficient (32), presently we need to maintain a proper balance 
between (small) noise intensity q and the (not too large) frequency uj c to guarantee a positivity of 
(69), thus allowing for a subsequent redefinition in terms of D 2 (t), c.f. Eqs. (32), (32) for comparison. 
If so, we can safely write (keep in mind that we consider exclusively the planar motion): 

'Aw 1 / 2 ' 



2(fg-h-k 2 )V 



,1/2 



-VQ 



(70) 



The important feature of this local conservation law is that on its right-hand side appears the 
very same force which is present in the Kramers equation l|49l) . as opposed to the corresponding 
conservation law in the Smoluchowski diffusion regime, |13j . 
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4 Harmonically bound particle 

4.1 Velocity-space process 

For completeness, let us analyze in detail the case of harmonically attracting force. We consider the 
one-dimensional harmonic oscillator with circular frequency u. The equation of motion reads (for 
further notational convenience we replace the previous A(t) by £(t)): 



dx 

— = u 
dt 

- = -utx+m 



(71) 



or 



(72) 



d_( x \ _ ( / 

dt\u)-\-u? 0){u) + { £(t) 

and differs from the Langevin equation in the standard Ornstein-Uhlenbeck theory in that we have 

no friction term —f3u. 

The formal solution to this equation is given in terms of the resolvent operator R (t, t ) — e A ^~ to \ 
1 
u? 



where A = ( 2 ~ I an d has the form 



~x (t) = R{t,t )~x + [ R(t,s)~i (s)ds . (73) 

For the sake of simplicity of the formulae we introduce the following notation ~x (t) = (x (t) ,u(t)) 
and in our case £ (s) = (0, £ (s)). 

The characteristic equation of the matrix A is A 2 + u 2 = 0, with the eigenvalues Ai = iu and 
A 2 = — iui. The corresponding eigenvectors are x 1 = ^ ^ ^ and x 2 = ^ ^ ^ . In the base of 

eigenvectors {a; 1 , a; 2 } matrix A has the diagonal form A' = ^ ^ 

4'f/-f / e Ai(*-*o) \ 

Since e l * °- 1 = I , e A 2 (t-t ) ) we can now easn y evaluate the resolvent operator 

R(t,t ) = e A( - t - t °'> = Be^^-^B- 1 , where B = ( 1 1 ] is the matrix related to the change 



Ai A2 

of base and B _1 = \ 2 ^L\ 1 ( ^\ 1 )' Namely, we have: 



1 / A 2 e Al <t-*o) - A ie A2 (*- t0 ) _ e Ai(t-to) +e A 2 (t-to) \ 

* ( ' oj ~ X^X; { A x A 2 e^(*- t o) - A 1 A 2 e A =(*-*o) -A ie M*-to) + A 2 e^<-*°) J (74j 



which can be explicitly written as follows: 



*(*,*>)=( cosw ^-^. £ sinw (*-*») , (75) 
! -wsinw (f — to) cosw(t — to) 
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The formal solution to (|71[) has the form 



1 r 1 

Xo cosw (t — to) + uq— sinw (t — to) + / — sinu (t — s) £ (s) ds 

— xoLusmuj (t — to) + uq cosw (t — to) + / cosw (t — s) £ (s) . 

"'to 



Let us introduce the following notation 



1 



R = x — xo coslu (t — to) — uo— siulu (t — to) 



S 



u + xqlu sin lu (t — to) — uq coslu (t — to) 



cp(s)£(s)ds 



(76) 
(77) 

(78) 
(79) 



to 



where ijj (s) = — sin a; (t — s) and 93 (s) = coscj (t — s). 
If we define furthermore: 



2q / V (s)rfs 



2? 



to 



sin u> (t — s) ds — — ^ (t — to) 

to w 
t ,-t 



2q tp (s) ds = 2q / cos lu (t — s) ds = q(t — to) 



1 



sin 2w (t — t ) 
2lu (t - 1 ) 

sin 2cj (t — to) 
2lu (t - t ) 



h = 2q tp (s) tp (s) ds = — / sinu (t — s) cos a; (t — s) ds = —r sin 2 ui (t — to) 

J to W J to w 



(80) 
(81) 
(82) 



we end up with a function W (i?, S) of the form (20) which can be read out as a transition probability 
density for the Markovian phase-space process (85), here denoted as P (x, v, t\ xo, vq, to), which clearly 
solves the appropriate Kramers-Fokker-Planck equation: 



dP dP 



dP d 2 P 



dt 



4.2 Local conservation laws 



-5 — I- v ~^~ = u x ~z, — I - c l 



dx 



dv dv 2 



(83) 



As an initial density we choose po(x,u) of Eq. (21) and propagate that according to p{x,v,t) = 
J P (x,u,t\x Q ,u ,t = 0) po (xo,v )dx du . 

Since P (x, u, t\ xq, uo, to) is the fundamental solution of the Kramers equation, the joint density 
p(x,u,t) is also the solution. This joint density can be put into the form p(x,u,t) = W (R,S), c.f. 
Eq. (96), where 



S = U + Xi n iUI Sin Ujt — Uini cos ojt 

R = x — Xi n i cos urt — Ui n i— smut 



ui(b 2 + 2qt + a 2 w 2 ) + (- (b 2 lu) + a 2 lu 3 ) cos(2tw) -q sin(2tw) 



2w 3 



lu (b 2 + 2 q t + a 2 lu 2 ) + lu (b 2 - a 2 lu 2 ) cos(2 1 lu) + q sin(2 1 lu) 



2lo 



(84) 
(85) 

(86) 

(87) 



sin(toj) (lu (b 2 — a 2 lu 2 ) cos(t lu) + q sin(t uj)) 
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The marginals p(x,t) = J p (x, u, t) du and p(u,t) — J p (x, u, t) dx read: 



p{x,t) 



2ir^ (u(b 2 + 2qt + a 2 w 2 ) + (- {b 2 w) + a 2 lj 3 ) cos(2 t w) - q sin(2 t lj)) 



x exp 



(x — Xi n i cos cut — Vi n i^ sin wtY 



2 (oj (b 2 + 2 q t + a 2 w 2 ) + (- (b 2 lu) + a 2 lj 3 ) cos(2 1 u) - q sin(2 1 w)) 



and 



p(u,t) 



\2irg) V \ 2. 9 



(w (fe 2 + 2qt + a 2 LU 2 ) + lj (b 2 - a 2 lj 2 ) cos(2tLu) + q sm(2tu>)) 



x exp 



(u + Xi n iUj sin Lot — Ui n i cos out) 



2±{lu {b 2 + 2qt + a 2 LU 2 ) + lu (b 2 - a 2 lj 2 ) cos(2iw) + q sin(2iw)) 



At this point we can evaluate < u(t) > 2 — [xi n iW sin(u>t) — Ui n i cos(ujt)] 2 and next: 



< u 2 (t) > - < u(t) > 



2 _ u{b 2 + 2qt + a 2 LU 2 ) + Lu(b 2 - a 2 L0 2 ) cos(2tuj) + qsm(2tuj) 

2uj 



(89) 



to have clearly exemplified the kinetic energy growth (e.g. < u 2 (t) > — < u(0) > 2 ). Notice that by 

taking the lj — > limit, we recover the force-free dynamics result (35). 

w 

J w(S,R)dS 



As before, let us introduce the auxiliary distribution W (S\R) = ,■ 'p^ c which has the standard 



Gaussian form 



W(S\R) = 



2tt 



We define local moments {u) x — J uWdu and (u 2 ) x = J u 2 Wdu so that 



-XiniLO Sin LOt + Uini COS + 



("\ = («>S+(ff-y) 



1 



^ — cos — Ui„i — sin o;t 



and local conservation laws follow: 



9w d ,. . . 

■ :«l £ I («>. 



-w 2 cc - 



dt x dx 

The local momentum conservation law Ij94(l takes the form 



<" 2 L - <«>' 



1 dw 

w dx 



(90) 

(91) 
(92) 

(93) 
(94) 



Aw 1 / 2 



,.1/2 



(95) 



18 



where the time-dependent coefficient reads: 



fg - h 2 



(96) 



-q 2 + 2b 2 qtu 2 + 2q 2 t 2 lo 2 + 2a 2 b 2 uj 4 + 2a 2 qtuj 4 + q 2 cos(2 tuj) + qu (-b 2 + a 2 lu 2 ) sin(2 1 lu) 



Again, a proper balance between the (small) intensity parameter q and the frequency lu must be 
maintained to secure the positivity of (96), since then only the notation D 2 (t) can be consistently 
introduced and the pressure-type contribution +VQ explicitly identified in Eq. (95). 

5 Kinetic theory viewpoint: role of "collision invariants" 

Previously, we have directly identified the microscopic motion scenario as the phase-space stochastic 
process to be followed by an individual particle. The usage of the hydrodynamical picture implicitly 
refers to the collective (ensemble) "gas" picture with its corresponding Boltzmann kinetic theory. 
Within the Boltzmann framework, Fokker-Planck-Kramers equations are derived under various sim- 
plifying assumptions, |12l l7). and effectively play the role of oversimplified kinetic equations, with the 
Boltzmann collision integral being replaced by suitable differentiable expressions. 

In particular, let us notice that all hither to considered Kramers equations have the general kinetic 
form: 



where / = f(x, u t t) is a sought for phase-space probability density while C(f) is a substitute for 
the Boltzmann collision integral, m being a mass parameter. 

Would we have worked with the " normal" kinetic equation in the Boltzmann form, then a standard 
method to produce local conservation laws is known to employ the collision invariants and heavily 
relies on an explicit form of C(f) which must yield: J C(f)d 3 u = 0, J ~uC(f)d 3 u = and in addition 
J ~u 2 C(f)d 3 u — 0. That implies respectively mass, momentum and energy conservation laws, as a 
consequence of the existence of microscopically conserved additive quantities, ^JEIi c -f- also the 
very recent paper \6'6\ . 

A concrete form of the "collision integral" C(f) while adopted to the considered before stochastic 
processes, can be easily deduced from the corresponding Fokker-Planck-Kramers equations. For exam- 
ple, the frictional dynamics involves C(f) = f?V^ + fi~u ■ V / while by omitting the /3-dependent 
term we are left with the frictionless variant C(f) — qV 2 -*/. 

It is important to notice that only J C(f)d 3 u = holds true identically for both frictional and 
frictionless cases. The frictional case does not respect two other kinetic identities. On the other hand, 
the frictionless dynamics respects one more identity, namely J ~uC(f)d 3 u — , but does violate the 
third one related to the energy conservation rule. 



2lo 




(97) 
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As an almost trivial remark we may consider an observation that one more step in the hierarchy of 
gradually more constrained kinetic problems may be completed by choosing suitable solutions of the 
Liouville equation, which we obtain by passing to the q — > limit in Eqs. (14), (37) and (83). That 
corresponds to replacing the stochastically forced dynamical systems by deterministic Newton laws, 
c.f. (13), (34), (71). Since for such solutions, the right-hand-side of Eq. (97) identically vanishes, we 
obviously have J ~u 2 C(f)d 3 u = 0. 

For those Liouville-type probability densities, all three collision invariants do vanish, satisfying the 
demands put forward in Ref. |33j to support the "idea that the classical kinetics could be underlying 
quantum mechanics" . The point is that there is manifestly no collisional kinetics involved in this 
particular case although all minimal requirements of Ref. |33j are fulfilled. 

By performing q — > limits of previously considered joint velocity-space probability densities we 
pass smoothly to deterministic solutions of the respective Liouville equations in the familiar form of 
the positive Wigner functions (up to notation adjustments, see e.g. |27||2H|), with rather obvious 
hydrodynamical connotations. 

As an illustrative example of this reasoning let us reproduce in detail the q — > limit of the 
frictionless dynamics of Section 2. We have (tilde is used to discriminate between q — and q ^ 
cases): 

= dpp^ ( JjL r - ""ST" ' ) ■ m 

Here, the normalization is preserved, since we have J (J p(x, u, t) dx) du = J (J p(x, u, t) duj dx — 
1, while for the marginals p(x, t) — J p(x, u, t) du and p (u, t) — J p(x, u, t) dx, in the very same q — > 
limit we get: 



and 



P(M)= U(a 2 + ^ 2 ) J ^ 2(a?+W) J" (1 ° 0) 



In order to find local moments (u) x — j u p ~^^ du and {u 2 ) x — J u 2 p ~ c ^^ du of the joint proba- 
bility density p(x, u, t) we compute 



p(x,u,t) _ /a 2 + b 2 t 2 \ 2 
p(x,t) ~ \ 2na 2 b 2 J 6XP 



V 



^U Ui n i a 2_j_^2^2 {x ^int ^ini^)^ 



o a 2 b 2 
Z a7+E 2 ~t 2 ~ 



(101) 



Consequently, the local moments are given by 



b 2 t 

(u) = U ini + (x - X ini - U ini t) (102) 

a 1 + b z t l 
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(u 2 ) x = (u)l + 



2 h 2 



a z h 



b 2 t 2 



and the force- free Euler-type equation in the form (8), c.f. is clearly recovered 



Of 



+ (u) V (u) =2a 2 6 2 V 



Ap(a;,t ) 1/2 



1/2 



= VQ 



(103) 



(104) 



where (c.f. eq. (10)) we may identify D 2 = a 2 b 2 with a squared diffusion constant, |10llllj . By setting 
formally D — we recover the standard quantum mechanical " hydrodynamics" , however with no 
intervention of any nontrivial kinetics, contrary to the hypothesis of Ref. |33| . 

Quite analogously, the planar motion in a uniform magnetic field gives rise to the Euler-type 
equation 

"Ap 1 / 2 " 



x B + 2a 2 b 2 V 



p 1 / 2 



F 



VQ 



(105) 



in conformity with the analogous result recovered in Eq. (69). Here, F denotes the Lorentz force 
contribution. The harmonic case implies: 



d d 



ApV2 



(106) 



again with no genuine kinetics involved. 



6 Concluding remarks 

We have investigated in detail exemplary solvable cases of the frictionlcss random dynamics. All of 
them share a common feature to yield the local momentum conservation law in the specific form 



F -> 
Tt)-* = — + 2d(t)V 



Aw 1 / 2 
w 1 / 2 



(107) 



where F denotes external force acting on the particle, and d = (det C) ~ where C is the covariance 
matrix of random variables (vectors) S and R (defined for each system) and n stands for the dimension 
of configuration space of appropriate system. 

1. free particle: R — x, F = 0, n = 1 

2. charged particle in a constant magnetic field: R = (x, y), F = ^ (^ > )tJ x B , n = 2 

3. harmonically bound particle: R — x, F — —mcu 2 x, n = 1 

In case of harmonic and magnetic confinement, we have identified parameter range regimes that 
allow for a positivity of the time dependent coefficient d(t) = D 2 (t) (in the force-free case it is positive 
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with no reservations), when the pressure-type contribution in Eq. (107) acquires a characteristic form 

of-^? = +VQ. 

By investigating the frictionlcss diffusion scenario we have given support to the validity of Eulcr- 
type momentum conservation laws in a carefully defined stochastic process context. In particular, 
that resolves problems with a consistent theoretical framework for diffusion processes which would 
show up Lorentz force effects in the local mean (e.g. rotational dynamics) for long times (exceeding 
a single period, c.f ^3122]). However, the frictionless motion needs developing methods of control for 
the otherwise untamed mean energy growth (energization phenomenon). 

As a byproduct of the discussion we have demonstrated that the usage of three traditional " collision 
invariants" representing microscopically conserved additive quantities does not produce a sufficiently 
general kinetic theory background for the derivation of the probabilistic counterpart of the Schrodinger 
picture quantum dynamics, at variance with the recent proposal of Ref. [22j- The previous discussion 
indicates as well that neither dissipative nor non-dissipative stochastic phase-space processes based 
on the white-noise kinetics are valid candidates to that end. 
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